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^ ! Abstract 

^ ■ We obtain the exact expression for the Von Neumann entropy for 

^ ■ an infinite bipartition of the XYZ model, by connecting its reduced 

qq | density matrix to the corner transfer matrix of the eight vertex model. 

CN ' Then we consider the anisotropic scaling limit of the XYZ chain that 

yields the 1+1 dimensional sine-Gordon model. We present the for- 
mula for the entanglement entropy of the latter, which has the struc- 
ture of a dominant logarithmic term plus a constant, in agreement 
with what is generally expected for a massive quantum field theory. 
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1 Introduction 



Recently, the investigation of entanglement properties in systems with many 
degrees of freedom has become a challenging line of research which has put 
into evidence a growing number of connections between quantum information 
and computational science, statistical mechanics, quantum field theory as 
well as formerly very far away topics, such as spin systems and black hole 
physics. 

In particular, much attention has been devoted to the study of strongly 
correlated spin and/or electronic models in one dimension, with the aim to 
unveil the relationship between entanglement properties of the ground state 
and quantum phase transitions. As soon as one has to deal with systems 
made up of more than two two-level subsystems, there does not exist a unique 
way of characterizing the degree of entanglement stored by the system itself. 
Thus, several different measures of entanglement have been proposed in lit- 
erature over the last few years and used to study the critical properties of the 
above mentioned systems. These studies include more standard definitions 
of entanglement measures, such as concurrence [U [3], as well as others, 
such as Renyi entropies [4], local [5] indicators, or fidelity |6j. 

In the realm of quantum field theory, one usually studies entanglement 
properties of a system via the computation of the so called Von Neumann 
entropy S = — Tr[p^ log Pa] related to the reduced density matrix pa of a 
subsystem A, as proposed in a seminal paper by Calabrese and Cardy [7], 
in which both the critical (conformal) and the free massive cases have been 
examined. As for the massless situation, much evidence has been accumu- 
lated proving that, for a wide class of lattice systems, from the formula for 
5* one can extract the central charge of its scaling limit conformal field the- 
ory. This is the case, for example, of spin k/2 XXZ-chains (k = 1,2,3...) 
[8] or of the SU(3) AFM Heisenberg chain [9]. In addition, some attention 
has been devoted to the relationship between entanglement properties of the 
vacuum and the boundary conditions of the theory [H HQ]. Important an- 
alytical progress has also been made in the context of massive integrable 
quantum field models, for which a framework for the computation of Von 
Neumann entropy has been developed using factorized scattering and form 
factor techniques [TT1 [T2]. 

In this paper we will first study the XYZ model and compute the exact 
expression for the entanglement entropy for an infinite bipartition of such 
quantum spin chain. As a measure of the entanglement, we will use the 
Von Neumann entropy of the density matrix associated with a semi-infinite 
segment of the chain. We will avoid the difficulties concerning the direct 
computation of the density matrix by mapping the quantum chain onto a 



2 



two-dimensional classical spin system. As Nishino et al. [13], [H] first pointed 
out, there is a connection between the density matrix of a quantum chain 
and the classical partition function of a two-dimensional strip with a cut 
perpendicular to it pj)]. In fact, the ground state of the quantum chain 
described by a Hamiltonian H is also an eigenvalue of the row-to-row transfer 
matrix T of the classical model, provided that [H, T] = 0. In [7j, this analogy 
was used in order to compute the entanglement of a transverse Ising chain 
and of a XXZ chain near their critical points. By exploiting the existing 
link between the eight vertex model and the XYZ chain [IB], we will obtain a 
formula for the entanglement entropy for the latter, which once more confirms 
the following universal expression for S, which is valid near criticality when 
the correlation length £ is much larger than the lattice spacing a: 

S ~ (c/6) log(£/a) + U (1) 

where c is the central charge of the conformal field theory describing the 
critical point that the theory is approaching. This formula first appeared in 
the context of black hole physics [T7] . 

Here U is a non-universal constant which is well known to depend on the 
particular model under investigation. When the bipartition is made up of two 
semi-infinite chains, as it is in this paper, it has been noted by several authors 
El EH HE] that it contains information about the so-called Affleck-Ludwig 
boundary entropy [19|. 

At the present time, with the method of this paper, we are not able to 
extract the boundary entropy information from such non-universal U. The 
study of the exact link between this term and the boundary entropy needs 
more accurate calculations on correlation functions that are not yet easily 
accessible in integrable models. The present exact result, however, could 
become very useful the day one will be able to compare it with calculations 
coming from independent methods. 

Finally, we will be interested in a particular and relevant scaling limit 
of the XYZ chain, which yields the 1+1 dimensional sine-Gordon model 
[20l[2T]. After a brief discussion of the connection between these two models 
in the thermodynamic limit, we will present the formula for the entanglement 
entropy of the latter which has a dominant logarithmic term in perfect agree- 
ment with what we expected by seeing the sine-Gordon model as a perturbed 
c = 1 conformal field theory. This formula also gives an analytic expression 
for the constant U. 
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2 Entanglement entropy for the XYZ model via 
corner transfer matrix 



Let us consider the quantum spin-| XYZ chain, which is described by the 
following hamiltonian 

H X YZ = -J2( J =° a n a n+l + Vn^l+l + J z a n a n+l) ( 2 ) 
n 

where the a l n (i = x, y, z) are Pauli matrices acting on the site n, the sum 
ranges over all sites n of the chain and the constants J x , J y and J z take 
into account the degree of anisotropy of the model. Without any loss of 
generality, we can put J x — 1. In the following we will exploit the very well 
known connection between the XYZ model and the eight vertex model [22]. 
Indeed, as shown by Sutherland [16], when the coupling constants of the 
XYZ model are related to the parameters T and A of the eight vertex model 
at zero applied field by the relations 

J x : J y : J z = 1 : T : A (3) 

the row-to-row transfer matrix T of the latter model commutes with the 
Hamiltonian H of the former. It is customary [22] to parametrize the con- 
stants A and Y in terms of elliptic functions 

l + fc 2 sn 2 (iA) cn(q)dn(zA) 

1 - k sn 2 (zA) ' 1 - k sn 2 (zA) 1 ' 

where sn(x), cn(x) and dn(x) are Jacobian elliptic functions while A and k 
are parameters whose domains are the following 

< A; < 1 , 0<A< I(k') (5) 

I(k') being the complete elliptic integral of the first kind of argument k! = 
Vl — k 2 . We recall that this parametrization is particularly suitable to de- 
scribe the anti-ferroelectric phase of the eight vertex model, that is when 
A < — 1, even if it can be used in all cases by redefining the relations that 
hold between T and A and the Boltzmann weights. 

Then, if relation ([3]) holds, in the thermodynamic limit one can obtain 
[T3l [T4l [T5] the reduced density matrix pi of the XYZ model relative to a 
semi-infinite chain as product of the four Corner Transfer Matrices A — C 

1 On the relation among T, A and the Boltzmann weights of XYZ we adopt here the 
conventions of [22J. 
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and B = D of the eight vertex model at zero field. The starting point is the 
density matrix defined by p —\ }( | , where | ) G H is the ground state in 
the infinite tensor product space H = 7i\ ® TL2 of our quantum spin chain. 
The subscripts 1 and 2 refer to the semi-infinite left and right chains, pi is 
then defined as the partial trace: 

Pi = Tr Wa (p) (6) 

More precisely one can write 

p 1 (a 1 a f ) = (ABCD)^, = (AB)l,, (7) 

where a and a' are particular spin configurations of the semi-infinte chain. 
Generally speaking, the above quantity does not satisfy the constraint Trpi = 
1, so we must consider a normalized version p[. Let us definine for future 
convenience the partition function Z = TrpV For the zero field eight ver- 
tex model with fixed boundary conditions, as it is shown in [22], [23J, EI], 
one can write down an explicit form of the Corner Transfer Matrices in the 
thermodynamic limit 

A{u) = C(u) = 

B(u) = D(u)=^ Q "J®^ " Js^ "J®... (8) 

where s and t are functions of the parameters A, k and u whose explicit 
expressions are given by s = exp [— iru/2I(k)] and t = exp [— 7r(A — u)/2I(k)), 
I(k) being an elliptic integral of the first kind of modulus k. Thus the density 
operator of formula ([71) is given by 




where x = (st) 2 = exp[— 7rA//(/c)]. We notice that p\ is a function of A and 
k only and that, furthermore, it can be rewritten as 

p 1 = (AB) 2 = e~ €d (10) 

O is an operator with integer eigenvalues and 

e = 7i\/I{k) (11) 

The Von Neumann entropy S can then be easily calculated according to 

S = -Trp[ In p[ = -e^rr^ + In Z (12) 
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where, in our case, the partition function is given by 

oo oo 

z= n(i+^) = n( i + e ^ W) ) ( i3 ) 

Thus we obtain an exact analytic expression for the entanglement entropy of 
the XYZ model 

oo oo 

Sxvz = e£ ^-y + £ Ml + e"^) (14) 

which is valid for generic values of A and k. This is the main result of 
this section. When e C 1, i.e. in the scaling limit analogous to the one 
of [8], formula (fT4l ) can be approximated by its Euler-Maclaurin asymptotic 
expansion, yielding 

7T 2 1 In 2 . 

= T7-- + 0(f) (15) 

This will be used in the next subsection where we will check our analytic 
results with some special known cases, namely the XXZ and the XY chain. 



3 Some checks against known results 

Let us first consider the case k = (i.e. r = 1) and the limit A — > + , which 
corresponds to the spin 1/2 XXZ. In this limit the eight vertex model reduces 
to the six vertex model. Let us note that formula f fT5l) coincides exactly with 
the one proposed by Weston [8] which was obtained in the study of more 
general spin k/2. In this limit the approximation e 1 is still valid and we 
can use the result of equation (TT5l) . From relation Q, it follows that 



A = v^V-A -1 + ((-1 - A) 3/2 ) (16) 
Thus equation (fT5l ) gives 

S = -^J— - ^ + O U-l - A) 1/2 ) (17) 

1 9,/9 -/-A _ 1 9 V 1 J y ' 



12V2 V-A - 1 

which can be written in a simpler form if we recall that, when e — > (i.e. 
A — > — 1 _ ), the correlation length is given by [22J 



c 2 

ln^ = — -21n(2) + 0(e~" 2/e ) 
a 
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Figure 1: Decoupled anisotropic Ising lattices. Horizontal and vertical lines 
belong to the original eight vertex model lattice, diagonal lines belong to the 
dual Ising lattice. 



where a is the lattice spacing. Recalling that 



e = 2V2V-A - 1 + O ((-1 - A) 3/2 ) 



(19) 



the expression for the entanglement entropy becomes 



S= Il n i + C/ + 0((-l- A) 1 / 2 



6 a 



(20) 



where U = — ln(2)/6. This last expression confirms the general theory of 
equation (pp) with c = 1, which is exactly what one should expect, the XXZ 
model along its critical line being a free massless bosonic field theory with 



As a second check, we consider the case T = 0, which corresponds to the 
XY chain. It is convenient now to describe the corresponding eight vertex 
model by using Ising-like variables which are located on the dual lattice [22J, 
thus obtaining an anisotropic Ising lattice, rotated by 7r/4 with respect to the 
original one, with interactions round the face with coupling constants J, J', 
as shown in figure [D 

In our case the Ising lattice decouples into two single sublattices with 
interactions among nearest neighbors. Now 



c = 1. 



A = sinh(2/?J) sinh(2/?J') = kj 



1 



(21) 



so that, using the elliptic parametrization, one has 



A = \l{k') 



(22) 




(23) 
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Let us approach the critical line of the anisotropic Ising model from the 
ferromagnetic phase, i.e. let us assume that kj — > 1~. In this case it is 
straightforward to write 



e = 



ln(l - fc/) 



+ O (ln- 2 (l - kj)) 



(24) 



so that the entanglement entropy is 



S = --\n(l-k I )+0(\n-\l-k I )) 



(25) 



Since £ 



(l-A; / ) + 0((l-fc / ) 2 



), we can easily conclude that 



S = Il n i + 0(ln- 1 (l-fc / 



6 a 



(26) 



where again the leading term confirms the general result (pp), with c = 1. 
This result is in agreement with what found in previous works (3j d [25] [26] 
by means of different approaches. 

4 The sine- Gordon limit 

In [20], ETj it has been shown that a particular scaling limit of the XYZ model 
yields the sine-Gordon theory. In this section we will use this fact to compute 
the exact entanglement entropy between two semi-infinite intervals of a 1+1 
dimensional sine-Gordon model. In his article, Luther [20] demonstrated 
that in the scaling limit, where a — > while keeping the mass gap constant, 
the parameters of the XYZ model and those of the sine-Gordon theory are 
connected by the following relation (keeping J x = 1 from the beginning) 




(27) 



where the parameter /j, is defined as 




(28) 



Here M is the sine-Gordon solitonic mass, and l r = I a ^ n , with 




(29) 
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These relations tell us how the coupling constant J z is connected to the 
parameter (3 of sine-Gordon, and how J y scales when we take the scaling 
limit a — > 0. It is clear from equation ( |29l ) that in this limit J y — > 1~. In the 
following we work in the repulsive regime 4n < (3 2 < 8n (which corresponds 
to0</i<7r/2 and — 1 < J z < 0). In this regime the mass gap of the theory 
is the soliton mass M. Taking this limit we use the following parametrization 
of the XYZ coupling constants 

T = y , A = | (30) 

Jy Jy 

which amounts to a reparametrization of the Boltzmann weights of XYZ suit- 
able for the | A | < 1 disordered regime where we are working now (see chapter 
10 of \'2'2\ for details). As a consequence of such reparametrization a minus 
sign appears in front of both equations (pj. Taking the sine-Gordon limit, A 
and k parametrizing V and A must now satisfy the following constraint 

Sn2(a) = - ( 31 ) 

Jy 

Considering the parametrization of A and using the properties of the Jaco- 
bian elliptic functions we can write 

_ cn 2 (.A)dn 2 (*A) _ - k% + (g - k% + % + l) 

(1 — k sn 2 (iA)) 2 ~~ (k 2 + k) 2 [ ' 

Expanding around k — ► 1~ and collecting A 2 = 1/J 2 from both sides of the 
equation we find 

A 2 = 1 + ^(1- J 2 ){k - if + 0{k - if (33) 
Using equation ( 1271) we obtain 



I 2 = iy*fr = A 2 ^ ( -^-) (34) 

\nsmnj 

where \i is completely fixed by choosing a particular value of J z . Now using 
the definition {2SD and fl3DD we find 

A 2 = 1 + (1 - J 2 )4 2 - 3 ^ ( W " + 0(a 4 ^) (35) 

Vvrsm/i/ v 1 
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which is valid when a — > 0. Comparing equation (1331) with (1351) we can 
identify in which way k scales to l - 

k = 1 _ 2 s(i-nM ( ^ + 0(a 2MA) (36 ) 

y 7r sin \i J 

Remembering the constraint ( 13TI ) and using the previous expression for k we 
have 

sn 2 (*A) = ~' Jz _~ l + O(a^) (37) 

When k — > 1 the elliptic function sn reduces to an hyperbolic tangent, thus 
we obtain 



tan 2 A = \——r + 0(a M/7r ) — ► A = arctan W j-^-y + 0(a^ n ) (38) 

Now we can evaluate the expression ( TT4l ) in this limit. Using the following 
asymptotic behaviour of the elliptic integral I(x) 

1 3 

I(x) w --ln(l-x) + -ln2 + 0(l-x), x « l - (39) 

along with the approximation ( fl~5l) . we can write the exact entanglement 
entropy of a bipartite XYZ model in the sine-Gordon limit 

7r ln(l — k) — 3 In 2 In 2 
Ss G = -WTT ~ — + °(Vln(a)) (40) 



arctan . 



1- J, 



The leading correction to this expression comes from the 0(e) term of equa- 
tion ( 1T51) . The constant J 2 is connected to j3 by 

J, = - COS 7^1-^ (41) 

thus using this property and the scaling expression Q3TH we can write down 
the entanglement entropy as 

1 / 1 \ 1 /sin 7r (1 - f^)l \ 

s «° = e ln Gss) + 5 ln ( (i-g) J + ln(0)) (42) 

This result confirms the general theory due to [ZlCEI], in the limit where the 
system is bipartite in two infinite intervals, with the central charge equal 
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to 1, as it should, because the sine-Gordon model can be considered a per- 
turbation of a c = 1 conformal field theory (described by a free massless 
boson compactified on a circle of radius \Ztt//3) by a relevant operator of 
(left) conformal dimension /3 2 /87r. We can write 



S, 



sG 



6 \Ma 



U{(5) 







(43) 



where the constant term U(f3) takes the value 



£/(/?) = -In 
o 



sin 


>(>-£) 











(44) 



At /? 2 = 47r, when the sine-Gordon model becomes the free Dirac fermion 
theory, it assumes the value C/(\/47r) = |log2 = 0.11552453..., while at 
(3 2 = 87r, where the theory becomes a relevant perturbation of the WZW 
conformal model of level 1 by its operator of left dimension |, it becomes 
C/(v/8tt) = | logvr = 0.19078814.... 

We notice that formula (l4*2i) yields the exact value of the overall constant 
U(/3), since it has been derived from equation (fl~5l) which is exact up to terms 
0(e). As mentioned in the introduction, U contains a contribution from the 
Affleck-Ludwig boundary entropy that we are not able to determine at this 
stage. It would be extremely interesting to be able to perform the same 
calculation in other regularization schemes of the sine-Gordon model. Also, 
consideration of boundary conditions other than the "vacuum" one consid- 
ered here could shed light, in the spirit of j8], on the problem of clearly 
relating this constant to the boundary entropy. In this respect the consider- 
ations done in paper [H] are very important. Clearly this question deserves 
further investigation that we plan for the future. 



5 Conclusions and outlooks 

We have carried out an exact formula for the entanglement entropy in an 
infinite bipartition of the XYZ model in the thermodynamic limit, by which, 
as test bench, we have re-obtained some well known results about the en- 
tanglement entropy of two integrable systems, the XXZ model and the XY 
chain. 

In addition we have obtained the entanglement entropy of the (repulsive) 
sine-Gordon scaling limit of the XYZ model, which, on one side, confirms 
once more the general theory due to [H [H] and, on the other side, yields 
for the first time an exact expression for such quantity in the case of an 
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interacting massive field theory. Since we are dealing with a non free theory, 
an investigation of the connections between this expression, Affleck-Ludwig 
boundary entropy [19] and one point functions of sine-Gordon fields [27] is 
non-trivial and deserves further work. 

Also it would be important to compare our results with those obtained 
for massive 1+1 dimensional theories by the form factor approach of [HI [12] 
and generalized to any massive (even not integrable) theory in [28]. This 
implies to implement our technique for a subsystem A consisting of a finite 
interval, a situation that we plan to study in the future. 

Finally, let us remark that another issue addressed in [7] is the dependence 
of entanglement entropy on finite size effects. To investigate this aspect one 
should implement finite size effects into the XYZ model, possibly within the 
Bethe ansatz approach, and then rescale to the continuum. A good question 
to ask is how far this approach can be related to the definition of sine-Gordon 
model on a cylinder by rescaling lattice models, both in the XYZ approach 
[29], [30] or in the so called light-cone approach by Destri and De Vega [31 J. 
How far the finite size effects on entanglement entropy may be encoded in 
structures generalizing the non-linear integral equations governing e.g. the 
finite size rescaling of the central charge c (see e.g. [32] and references therein) 
and the Affleck-Ludwig boundary entropy g [33, 34|, is an intriguing issue, 
surely deserving further investigation. 

All these lines of developments can shed new light on the comprehension 
of the entanglement problem in quantum field theory. 
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